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ABSTRACT 

We examine matter-enhanced neutrino flavor transformation {i^t{^i) ^ ^e) in the re- 
gion above the neutrino sphere in Type II supernovae. Our treatment explicitly includes 
contributions to the neutrino-propagation Hamiltonian from neutrino-neutrino forward 
scattering. A proper inclusion of these contributions shows that they have a completely 
negligible effect on the range of i^e-^ri^j.) vacuum mass-squared difference, 5m^, and vacuum 
mixing angle, 6, or equivalently sin^ 26, required for enhanced supernova shock re-heating. 
When neutrino background effects are included, we find that r-process nucleosynthesis 
from neutrino-heated supernova ejecta remains a sensitive probe of the mixing between a 
light Uf, and a z^r(/i) with a cosmologically significant mass. Neutrino- neutrino scattering 
contributions are found to have a generally small effect on the (Jm^, sin^ 2^) parameter 
region probed by r-process nucleosynthesis. We point out that the nonlinear effects of the 
neutrino background extend the range of sensitivity of r-process nucleosynthesis to smaller 
values of 6m'^. 
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I. Introduction 



In this paper we investigate the problem of matter-enhanced neutrino flavor transfor- 
mation in the region above the neutrino sphere in Type II supernovae. In particular, we 
examine the role of contributions to the neutrino-propagation Hamiltonian from neutrino- 
neutrino forward scattering. A general framework for treating these contributions in the 
context of the Mikeheyev-Smirnov-Wolfenstein (MSW) neutrino flavor transformation pro- 
cess has been given in Ref. [1] (see Ref. [2] for a numerical study of the case of a pure 
neutrino gas). Although the role of neutrino- neutrino scattering in the problem of matter- 
enhanced neutrino flavor conversion in supernovae has been treated previously [3, 4], the 
present paper gives the flrst complete treatment utilizing the scheme of Ref. [1] . 

Recent studies have examined MSW tranformation of Vr or into Ve in the region 
above the neutrino sphere in the post-core-bounce supernova environment [5, 6]. These 
studies suggest that if v^. or has a mass in the cosmologically interesting range of 1-100 
eV, then matter-enhanced transformation to will be possible in this region. Such trans- 
formation can result in significant effects on supernova dynamics and/or nucleosynthesis. 

If we define, for example, \ve) and \vr) to be flavor eigenstates of Vf. and Vri and 
and \v2) to be the associated mass eigenstates, then the vacuum mixing angle, 9, is deflned 
through 

l^-e) = cos ^jz/i) -h sin ^11/2), (lo) 

= —sm.9\vi) -\- cos9\v2). (16) 
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Ref. [5] shows that i^t{h) ^ mixing with sin^ 29 > 10~^ in the region above the 
neutrino sphere at a few hundred milhseconds after the bounce of the core can result in a 
(30-60)% increase in the supernova shock energy. Ref. [6] shows that the heavy element 
nucleosynthesis from the hot bubble region is sensitive to i^t{ij,) ^ mixing at a level of 
sin^ 29 ~ 10~^. This hot bubble region forms above the neutrino sphere ~ 3 seconds after 
core bounce. These effects are sensitive to mixing angles far smaller than those which can 
be probed in laboratory experiments. These supernova effects ultimately may represent 
our most sensitive probe of putative neutrino dark matter. 

However, studies [5] and [6] neglected the off-diagonal contributions of neutrino- 
neutrino scattering to the flavor-basis neutrino-propagation Hamiltonian. In what follows, 
we present a detailed study of neutrino flavor transformation in the post-core-bounce su- 
pernova environment. Our calculations include all effects of the neutrino background. We 
have adopted the overall principles and techniques of Ref. [1] in our treatment of neutrino- 
neutrino and neutrino-electron scattering contributions to the neutrino-propagation Hamil- 
tonian. We flnd that neutrino background contributions have a negligible effect on the 
range of Ve-^riix) vacuum mass-squared difference, 5m^, and vacuum mixing angle, 9 (or 
sin^2^), required for enhanced supernova shock re- heating. A proper treatment of the 
ensemble average over the neutrino background shows that r-process nucleosynthesis from 
neutrino-heated supernova ejecta remains a sensitive probe of the mixing between a light 
z^e and a v^- (or v^j) with a cosmologically signiflcant mass {jn^^^^^ ~ 1-100 eV). 

In Sec. II we discuss a general framework for treating neutrino flavor transformation 
in the supernova environment. In Sec. Ill we compute neutrino flavor transformation 
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probabilities as functions of 5m^ and sin^ 26 relevant for the shock re-heating and hot 
bubble/r-process nucleosynthesis epochs of the supernova. We give conclusions in Sect. 
IV. 

II. The Neutrino-Propagation Hamiltonian in Supernovae 

The general problem of the time evolution of the full density matrix for an ensemble 
of three flavors of neutrinos and antineutrinos with an electron/positron background and 
a nucleon background is a daunting one. Several formal approaches to this problem have 
been made (cf. Ref. [1] and references therein). In the present paper, we shall only 
summarize the salient features of this previous work and taylor our subsequent discussion 
to the particular problem of neutrino propagation and flavor transformation in the region 
of the supernova environment above the neutrino sphere. Considerable simplification of 
the problem can be realized in this case. 

The general time evolution of the neutrino density matrix p can be summarized as 

ip=[H,p], (2) 

where p = Pij\'i){j\, P = dp/dt, and i and j refer to all neutrino quantum numbers 
including momentum (energy), flavor, helicity, charge conjugation eigenvalue, etc. In Eq. 
(2), H is the full Schrodinger picture Hamiltonian including all neutrino self interactions 
as well as interactions with the and nucleon backgrounds. 

Without loss of generality we can follow a particular momentum component of Eq. (2) 
(cf. Ref. [1]), or equivalently, the associated Schrodinger equation for the time evolution 
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of neutrino field ampfitudes for a given momentum. The Hamiltonian operator in this case 
would have the dimensionality of the density matrix for the single momentum state (e.g., 
12x12 for three Dirac neutrino flavors, since each neutrino state has either right-handed 
or left-handed helicity, and is either a neutrino or an antineutrino). 

We argue that further simplification of this problem can be made through approxima- 
tions motivated by the par^zca/ar distribution functions for z^e: ^e-, ^/j., ^n-, and Vt which 
obtain in the region above the neutrino sphere in the post-core-bounce epoch of Type II 
supernovae. Since the distribution functions for zv^, iv^, and V-j- are all expected to be 
essentially identical, mixings between neutrinos in this sector will have no effect on any 
aspect of supernova physics. In other words, we need only consider mixings between v^. 
and either or v^-. If, as seems likely, the vacuum mass heirachy for neutrinos satisfies 
"^''t(m) ^ '^v^-i then we need only consider matter-enhanced mixing among neutrinos, as 
antineutrino mixing is supressed by matter effects. 

The masses "t-i^^^^^ 1-100 eV of interest in the post-core-bounce supernova environ- 
ment are very small compared to the typical neutrino energies (average neutrino energy 
{Ejy) is about or greater than 10 MeV). In this case we can neglect the population of 
right-handed Dirac neutrinos and left-handed Dirac antineutrinos produced by scattering 
processes. This is because helicity-flipping rates are proportional to {rriy/EyY. 

Taking advantage of these features allows us to reduce the dimensionality of the Hamil- 
tonian in Eq. (2) to 2 x 2 for the Dirac neutrino case. If neutrinos are Majorana particles, 
then we have only left-handed neutrinos and right-handed antineutrinos, and again the 
Hamiltonian of interest is 2 x 2. 
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In any case, the neutrinos of interest in supernovae will be extremely relativistic, so 
that we can approximate the neutrino energy as Ei^ = \J'p^ + ~ p + w? /2p. The first 
term in this expression, p, the momentum, just gives an overall phase to the coherent 
propagating neutrino state and can be ignored without loss of generality. The second 
term, m? /2p, is responsible for the relevant neutrino mixing behavior. The part of the 
Hamiltonian corresponding to the /2p term in vacuum, H^, can be written in the flavor 
basis (e.g., \ve)i l^r)) as 

^ ~ 2 V sin2^ cos2^ ) ' 
where 9 is the vacuum mixing angle as in Eq. (1), A = 5w? /2Ey, and 5m^ = m| — m\, 
with 777-1 and m2 the vacuum mass eigenvalues corresponding to the mass eigenstates \vi) 
and \i'2)i respectively. 

In matter the relation between the flavor basis and the mass basis can be written as 
in Eq. (1), but with the vacuum mixing angle replaced by an appropriate matter mixing 
angle On- For illustrative purposes consider the case where the only contribution to the 
effective mass difference between neutrino flavors comes from charged-current exchange 
scattering on electrons. We take the net number density of electrons to be 

rte = ng- — ne+ , (4a) 

where n^- (ne+) is the total proper number density of negatrons (positrons). The electron 
fraction Yg. is deflned in terms of the total baryon rest mass density p and Avogardro's 
number Np^ by 

Tie = pYeNA. (46) 
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The contribution to the Hamiltonian from neutrino-electron exchange scattering is 

A = \/2G'Fne, (5) 

where Gp is the Fermi constant. In Fig. 1 we show a generic Feynman graph for v^,- 
e scattering. To obtain the result in Eq. (5) one must sum graphs for i^e-e~ and fe- 
e'^ scattering over the appropriate distribution functions. In this case the neutrino- 
propagation Hamiltonian, i/g) can be written as 

_Aeff/-cos2^„ sin20^\ _ 1 /-Acos20 + A Asin2^ \ 

2 V sin26'n cos2^„j~2V Asin26' Acos26'-Aj' ^> 



where Aefi = y ( A cos 29 — + A^ sin^ 29. In these expressions the matter mixing angle, 
9ni is related to the vacuum mixing angle 9 and the local net electron number density 
through 

. c^r. Asin2^ 

sm 29n = I = , (7a) 

^(A cos 29-Ay + A^ sin^ 29 

Acos29-A 

cos29n = ■ (76) 

^(A cos 26* - A)2 + A2 sin^ 29 

The amplitudes for antineutrino-electron (Pg-e) exchange scattering and neutrino- 
electron {ue-e) exchange scattering have opposite signs. This implies that Ug-e exchange 
scattering gives a contribution —A to the flavor-basis interaction Hamiltonian for Pg- In 
this case the matter mixing angle for antineutrinos, 9n, satisfies 

. Asin26' ,^ , 

sm 29n = , (8a) 

^J (A cos 2^ + yl)2 + A2 sin^ 29 
a cos 26* a 

cos2^„ = . (86) 



y/(Acos2^T^4)2TA2^h^^ 



We note that the vacuum mixing angles for the neutrino and antineutrino sectors are the 
same. 

It is evident from Eqs. (7a & b) that matter effects can give enhancement of flavor 
mixing in the neutrino sector. Mixing is maximal at a mass level crossing, or resonance, 
where A cos 2^ = A [7]. On the other hand, Eqs. (8a & b) show that matter effects give 
supression of flavor mixing in the antineutrino sector. 

In the supernova environment, however, the neutrino background and the resultant 
neutrino-neutrino forward exchange-scattering effects necessitate some modification of the 
above treatment of neutrino flavor transformation. In the region above the neutrino sphere 
in post-core-bounce Type II supernovae the neutrino fluxes can be sizable (see, for exam- 
ple, the discussion in Ref. [6]). Individual neutrinos emitted from the neutrino sphere 
can be described as coherent states. However, each emitted neutrino is related to every 
other emitted neutrino in an incoherent fashion. In other words, these different individual 
(or single) neutrino states have random relative phases, as is characteristic of a thermal 
emission process. The total neutrino field is properly a mixed ensemble of individual neu- 
trino states. It is not a coherent many-body state. Accordingly, the total neutrino density 
matrix is an incoherent sum over each single neutrino density matrix. 

For a single neutrino emitted at the neutrino sphere as a Va (e.g., in fiavor state 
a = e, T for the case of two-neutrino mixing) we can represent its state at some point 
above the neutrino sphere as 

= aiait)\ui{t)) + a2ait)\u2it)), (9) 
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where li'iit)) and \i'2it)) are the instantaneous physical mass eigenstates of the full 
neutrino- propagation Hamiltonian, and aia{t) and a2a(^) are the corresponding complex 
amplitudes. Normalization requires that we take |aia(t)P + |a2a(^)P = 1. In these expres- 
sions the time, t, could be any evolutionary parameter (e.g., density, radius, etc.) along 
the neutrino's path from its creation position at the neutrino sphere to a point at radius 
r. The single neutrino density matrix is then given by 

|V'.J(V'.J = \aia{t)\yi{t)){Mt)\ + \a2ait)\y2{t)){Mt)\ 

+ aUt)a2a{t)\u2{t)){u^{t)\ + a,a{t)aUt)\Mt))Mt)\- (10) 

The density matrix representing the mixed ensemble of single neutrino states all with 
momentum p can be written as the incoherent sum 

Ppd^P = ^dn,,JiP,,J{iP,,J. (11) 

a 

In this expression the sum runs over, for example, a = e, t, while dn,y^ is the local 
differential number density of neutrinos with momentum p in interval d^p. The local 
diflFerential Va neutrino number density at a point at radius r above a neutrino sphere with 
radius Ri, is 

dn.^ ^ nlj^SEuJdE^^ , (12a) 

where dO,^ is the differential solid angle (pencil of directions) along the neutrino momen- 
tum p (|p| fa -Ei/^), is the neutrino number density at the neutrino sphere, and 
fy^ {Ey^ ) is the normalized fa. energy distribution function. We can show [6] that a good 
approximation for is 
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where L^,^ is the luminosity in neutrinos, {E,^^) is the average Va neutrino energy, and 
c is the speed of hght. The normahzed neutrino energy distribution function can be 
well approximated by 



where the rank 2 Fermi integral with argument zero is ^2(0) ~ 1.803, and where T^,^ is 
the neutrino sphere temperature. The average neutrino energy is related to the 
appropriate neutrino sphere temperature by {E^J) w 2t.lbT^^. 

In the region of the supernova above the neutrino sphere, the range of the solid angle 
contribution allowed in Eq. (12a) is restricted to be within the solid angle subtended by 
the neutrino sphere as seen from a point at radius r. The geometrical arrangement of a 
neutrino sphere with radius a point above the neutrino sphere at radius r, and various 
neutrino paths are depicted in Fig. 2. 

We can now write the full flavor-basis neutrino-propagation Hamiltonian as a sum 
of vacuum mass and electron background contributions, i/g, and neutrino background 
contributions, H^,,^: 



where H^,,^ represents the ensemble average over neutrino-neutrino interactions using the 
density matrix in Eq. (11). For a neutrino with energy Ej^ and momentum p propagating 
radially outside the neutrino sphere we can write 




F2(0)T3^exp(i?.„/T.J + l' 



(12c) 



H = H, + H, 



(13a) 
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where pq is the density matrix for antineutrinos with momentum q (defined in obvious 
analogy to pq in Eq. [11]) and 9q is the angle between the direction of the propagating 
neutrino with momentum p and the directions of other neutrinos in the ensemble with 
momentum q. We can generalize the expression for Hi,i, in Eq. (13b) for non-radially 
propagating neutrinos by replacing cos6'q with q-p/{|q||p|}- 
It is convenient to recast Eq. (13b) in the form 

I f B B, 



2 (sl -S) + ^^'^ - ^°^^)TKPq - P^)d\ (14) 



Note that the second term in this equation is simply proportional to the identity matrix, 
implying that it provides only an overall phase in the propagating neutrino state and can 
be ignored. 

In the first term in Eq. (14) there are two contributions to the neutrino-propagation 
Hamiltonian, B and B^j. {B^-g), where 

B = V2Gf j {1- cos 6'q){(pq - Pq)ee - (Pq - Pq)rT}c^^q, (15a) 
Ber = 2V2GF y"(l - COS ^q) (pq - Pq)erC?^q, (156) 
Bre = 2V2GF J {I- COS^q)(pq - Pq).eC^^q, (15c) 

where, for example, by (pq)er we mean the matrix element of the density matrix operator, 
{uelPdli^r), while by (pq)er wc mean (zyelpql^'r)- 

Here B corresponds to the forward neutrino-neutrino exchange-scattering contribu- 
tions to the neutrino effective mass. These contributions are the analogs of the Ue-e 
exchange-scattering term. A, in Eqs. (5) and (6). Generic Feynman graphs for these 
neutrino-neutrino exchange processes are shown in Fig. 3a for Vg-J^e scattering and in Fig. 
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3b for fr-J^r scattering. We will later refer to B as the "diagonal" contribution of the 
neutrino background to the flavor-basis neutrino-propagation Hamiltonian. 

The neutrino background also provides "ofF-diagonal" terms in the flavor-basis 
neutrino-propagation Hamiltonian. These are, for example, the Ber and B^-e terms above. 
They arise because the background neutrinos are not in flavor eigenstates [1]. We show 
graphically these contributions for ly^ and i^j- neutrinos with momenta p and q in Fig. 4. 
The corresponding diagonal and off-diagonal contributions to the flavor-basis antineutrino- 
propagation Hamiltonian from the neutrino background are —B and —B^t {—Bre), respec- 
tively. 

Considerable simpliflcation in the evaluation of B, B^j. {B^-^) can be realized by adroit 
attention to the phases in the expression for the single neutrino density matrix in Eq. (10). 
Note that the last two terms in Eq. (10) are cross terms. They have coefficients al{t)a2{t) 
and ai{t)a\{t), respectively. Each cross term is proportional to a factor ~ exp[z J uj\2{t)dt\^ 
with a;i2 the difference in the local neutrino flavor-oscillation frequencies of the two mass 
eigenstates \i^i{t)) and \v2{t)). These oscillation frequencies are, in turn, dependent on the 
local density. 

In both the early post-core-bounce shock re-heating epoch (time post-core-bounce 
tpB ~ 0.1-1 s) and in the hot bubble/r-process nucleosynthesis epoch (tpB ~ 3-15 s) the 
electron number density predominantly determines the neutrino flavor-oscillation frequency 
in the region just above the neutrino sphere [5, 6]. This is because the net neutrino number 
densities are negligible compared to the electron number densities (tt-e ~ 10^^ cm~^) in 
this region [5, 6]. The Hamiltonian in Eq. (6) is by itself sufficient to determine the 
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neutrino flavor-oscillation frequencies in this region. The local neutrino flavor-oscillation 
frequency difference in this case is u)i2 ~ Ags = (A cos 26 — A)'^ + sin^ 26. 

Furthermore, for the cosmologically interesting range Sm'^ = 1-100 eV^, the electron 
number density near the neutrino sphere greatly exceeds the MSW resonance densities 
for neutrinos with energies Ei, ~ 10 MeV. In this case the neutrino paths shown in Fig. 
2 will always cross a region where the electron number densities dominate the neutrino 
number densities prior to entering the resonance region. Therefore the local neutrino 
flavor-oscillation frequency difference will be a;i2 ~ Aefj ~ ^ = \/2Gpne along some part 
of every neutrino's path. 

When taking the ensemble average over the neutrino background we necessarily in- 
tegrate pp over neutrino momentum directions to a point at radius r. We thereby also 
average over the oscillating cross terms in Eq. (10). In addition, neutrinos with different 
momentum directions travel on paths with different lengths to arrive at a point at radius 
r. These different path lengths then give rise to different phases for the oscillating factor 
~ exp[i J U!i2{t)dt] in the cross terms in Eq. (10). In fact it is clear from Fig. 2 that each 
neutrino path from the neutrino sphere to a point at radius r will have a path length which 
depends on the polar angle. For neutrinos with momentum magnitude |p| each path with 
a different polar angle will have a different phase entering into the cross term coefficients 
ofEq. (10). 

The phase difference d(p acquired in going through a region of electron number density 
Ue with a difference in path length Sr is then 

64> ^ V2GFneSr ^ 642 ( J" ) ( • (16) 

\10'^° cm~"^/ \1 cm J 
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Path length differences of order 1 cm give rise to phase differences of » 27r ! It is obvious 
that the cross terms in the single neutrino density matrix in Eq. (10) vanish when aver- 
aged over all neutrino momentum directions. This argument can also be applied to the 
evaluation of pp since lvi^ > V^GpUe everywhere for antineutrinos. 

Clearly, we need only consider the first two terms of Eq. (10) in evaluating matrix 
elements of the density matrices pp and pp (cf. Eq. [11]). This will allow considerable 
simplification in computation of B and (Sre) from Eqs. (15a-c). 

Failure to properly perform the angular part of the ensemble average would result in 
the retention of non-zero cross terms in the neutrino density matrix elements, Eqs. (15a- 
c). This would introduce a spurious, and unphysical, "coherence" in the treatment of the 
neutrino background. In fact, the angular part of the ensemble average over the neutrino 
background is a key point in determining neutrino fiavor evolution in the region above the 
neutrino sphere in supernovae. 

Note that = 3^-^ since the terms in the ensemble averages, Eqs. (15a-c), are 
all real and the Hamiltonian must be Hermitian. The terms in Eqs. (15a~c) are all real 
because of the vanishing of the cross terms in the momentum direction average over the 
single neutrino density matrices. The full fiavor-basis Hamiltonian which includes both 
the electron and neutrino backgrounds is now 

rr rr ^rr 1 f - A COS 29 + A + B A Siu 2^ + \ ..^x 

H = H, + H,,= -\^ Asin2^ + Se. Acos2e -A-b)' ^^^""^ 

In analogy to the discussion preceding Eq. (6) we can rewrite this Hamiltonian as 

Ah/-cos26'h sin26'H^ 
2 V sin2^H cos2^Hy " ^ ^ 
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In this expression we have defined a full effective mixing angle, On, which, in analogy to 
Eq. (1), gives the relations between the flavor basis and the instantaneous mass basis 
including the effects of both the electron and neutrino backgrounds: 

\ue{t)) = cos^H(t)ki(t)) +sin^H(t)k2(t)), (18a) 
\ur{t)) = - sin eii{t)\iyi{t))+ COS euit)\u2{t)). (186) 

We have deflned Ah as 



Ah = V(A cos 29- A- By + (A sin 29 + B^rf. (19) 

The full effective mixing angle satisfies 

. Asin2gH + ^er ... ^ 

sm2^H = — , (20a) 

v/(A cos 29-A-BY + (A sin 29 + B^rf 

Acos2^H-A-S 

cos2^H = ■ (206) 

V(A cos 29- A- BY + (A sin 29 + B^rY 

Note that in the absence of a neutrino background Ah = Aeff and 9u = 9n. The corre- 
sponding expressions for the full effective mixing angle, ^h, in the antineutrino sector are 
obtained by replacing A, B, and B^j. with —A, —B, and — Sgr, respectively. 

Since the cross terms in the single neutrino density matrix will give no contribution 
to the ensemble average, we can write a reduced expression for the single neutrino density 
matrix in terms of flavor-basis eigenbras and eigenkets: 

1 1 

i\'ipiyj{fpuj)reduced = ^2 ~ ^2 ~ \(^ia{t)\^] COs29ii{t)}\l^e) {l^e\ 
+ - + l«la(^)l'] COs29Kit)}\l^r){l^r\ 

+ - |aia(t)P]sin2^H(t)(ke)K| + kr)^e|). (21) 
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With this form for the single neutrino density matrix, it is straightforward to evaluate 
flavor-basis matrix elements of the density matrix operator. For example, the expressions 
in Eqs. (15a-c) become 

B = -V2GfY, f (l-cos^q){[l-2|ai«(t)|2]cos2^H(t)cin,„ 

- [1 - 2|aia(t)H cos2eYi{t)dnyJ, (22a) 

= V^GfJ] / (I - cos^q){[l " 2|ai«(t)P] sin2en{t)dn^^ 

- [1 - 2|ai«(t)n sui2eYi{t)dn^J. {22b) 

In these expressions aiait) is the amplitude to be in the instantaneous mass eigenstate 
\j^i{t)) for an individual neutrino of momentum q which was created at the neutrino sphere 
{t — 0) in flavor eigenstate |j^o,). Likewise, dia{t) is the amplitude to be in the instantaneous 
mass eigenstate for an antineutrino of momentum q created at the neutrino sphere 

in flavor eigenstate li^a). The expressions dn^^ and duj^^ are as given in Eq. (12a), e.g., 
drir.^ <J^AE.JdE,JdnjA7i). 

It remains to evaluate these expressions for the particular conditions (electron density 
run and neutrino distribution functions) which obtain for the shock re-heating and hot 
bubble/r-process nucleosynthesis epochs. 

III. Neutrino Flavor Transformation in the Supernova Environment 



In this section we examine neutrino flavor transformation in the region above the neu- 
trino sphere in models of post-core-bounce Type II supernovae. There are several aspects 
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of the problem of neutrino flavor transformation in supernovae which are significantly dif- 
ferent from conventional computations of MSW flavor conversion in the sun. Foremost 
among these is the necessity of treating the neutrino background. In addition, the geome- 
try of neutrino emission from a neutrino sphere in a supernova is quite different from the 
solar case, where the neutrino source is distributed throughout the core. 

Bearing these points in mind, we can formally transform the full flavor-basis Hamil- 
tonian in Eqs. (17a & b) to the basis of the instantaneous mass eigenstates \i^i{t)) and 
\i'2{t))- The Schrodinger equation for the time evolution of the amplitudes aiQ,(t) and 
(J'2a{t) (see Eq. [9]) in this basis is then, 

. f ai^{t) \ _ f-An{t)/2 -iO^it) \ ( ai«(t)\ , . 

V «2a {t))-\ iOn (t) Ah {t)/2 J \ a^a (t) J ' ^ 

where aia{t) = daia{t)/dt, a2ait) = da2ait)/dt, and ^h(^) = d$ii{t)/dt. In this expression 
we follow the treatment of neutrino propagation and flavor transformation in Ref. [8] . Eq. 
(23) represents a set of nonlinear Gist order differential equations for the amplitudes aice{t) 
and a2ait)- The nonlinearity arises since, in general. Ah and the full effective mixing angle 
^H each depend on the neutrino background contributions B and (Eqs. [19], [20a & 
b]). In turn, B and Ser depend on the amplitudes aice{t) as in Eqs. (22a & b). 

The time evolution of the full effective mixing angle can be found from Eqs. (20a Sz 
b) to be 

• ^ Ber{Acos2e-A-B) + {Asm2e + Ber)iA + B) 
2[{Acos2e-A-B)^ + {Asm2e + Ber)^] ' 

where A = dA/dt, B = dB/dt, and B^t = dBer/dt. 
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We can define an "adiabaticity parameter" ^{t): 

7(0 - (25) 
2|9H(t)l 

Clearly, the neutrino mass eigenstate evolution is well approximated as being adiabatic 
when 7(t) » 1. Of course, if = 0, the neutrino mass eigenstate evolution is completely 
adiabatic, as can be seen directly from Eq. (23). 

The adiabaticity parameter generally satisfies ^{t) » 1 well away from resonance 
regions (neutrino mass-level-crossing regions). However, neutrino flavor conversion proba- 
bilities depend crucially on ^{t) at resonance. We shall denote the value of the adiabaticity 
parameter at resonance as 7(tres)- Resonance occurs when 

Acos26» = ^ + S. (26) 

We denote the position of this level-crossing point, or resonance, by tres- At resonance, 

^ _ (Asin2g + BeO^ _ {Asm2e + Berf ^ d\n{A + B) 
\A + B\ ~ A cos 2^ I dt 

The Landau-Zener probability for the neutrino to jump from one mass eigenstate to 
the other in the course of transversing a resonance region is [8], 

Plz ~ exp[-|7(tres)]- (28) 

Unlike the case for solar neutrinos, this expression is always sufficient for calculating neu- 
trino flavor transformation in supernovae [5, 6]. The Landau-Zener formula Eq. (28) is 
inapplicable for solar neutrino flavor conversion when, for example, neutrinos are created 
close to their resonance positions. This never occurs in supernovae, where neutrinos orig- 
inate on the neutrino sphere. The neutrino sphere is always well away from the resonance 
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region for the cases we will consider. In addition, solar neutrinos can experience double 
level crossings when they are created at densities below their resonance density. This does 
not occur in the post-core-bounce supernova environment. 

The very small vacuum mixing angles we shall consider for neutrino flavor conversion 
in supernovae imply narrow resonance regions. Narrow resonance regions, together with 
the generally large density scale heights (0.5-50 km) characterisitic of the region above the 
neutrino sphere [5, 6], imply that the first order Landau-Zener jump probability expression 
in Eq. (28) is always adequate [5, 6]. By first order jump probability we mean that we 
approximate the density profile as linear across the resonance region. 

It is obvious in Eqs. (24-27) that we recover the pure electron-driven neutrino flavor 
conversion case when the neutrino background disappears (i.e., B and B^j. vanish every- 
where). The neutrino background influences neutrino flavor evolution through resonances 
in two ways. 

First, the diagonal contribution of the neutrino background, B, essentially shifts the 
position of the resonance from the case where only the electron contribution. A, is present. 
This is evident from Eq. (26). The diagonal contribution of the neutrino background also 
alters the density scale height of weak interaction scattering targets at resonance. The 
density scale height of weakly interacting targets {\dlnn/dr\~^ following Eq. [7] in Ref. 
[6]) is the |cZln(^ + B)/dt\-'^ term in Eq. (27). 

The off-diagonal contribution of the neutrino background, B^t, has the effect of al- 
tering the adiabaticity of the neutrino flavor evolution at resonance. This is clear from 
Eq. (27), where B er appears in the expression for 7 (ires)- If A sin 2^ ^ |-Ser| then the 
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ofF-diagonal neutrino background contribution will have little influence on the adiabaticity 
of neutrino flavor evolution. 

However, the diagonal and off-diagonal contributions of the neutrino background in- 
fluence neutrino flavor evolution in a nonlinear manner, as outlined above. Not only are 
B and Sgr determined by the local neutrino distribution functions, but the local neutrino 
distribution functions are also dependent, in general, on the detailed history of neutrino 
flavor transformation. 

The crux of the problem of treating the nonlinear effects of the neutrino background is 
the computation of B and for the particular local neutrino distribution functions which 
obtain in the supernova environment. This will be evident if we discuss a simple iterative 
procedure for computing neutrino flavor transformation at resonance in the presence of a 
neutrino background. 

We can employ the Landau-Zener transformation probability in Eq. (28) to estimate 
the neutrino flavor conversion probability for a neutrino propagating through a resonance 
with the following simple procedure. We choose a vacuum mass-squared difference Sm^ 
and a vacuum mixing angle 9 (equivalently, sin^ 29) for a propagating neutrino of energy 

(1.) To begin with, we assume that Ser = 0. We use Sm^ and sin^ 2^, along with 
Ber = 0, in Eqs. (20a & b) to get a zero-order estimate for cos2^h, sin2^H, cos2^h, 
and sin2^H- Note that the value of A and B which enter into the expressions for cos26'h, 
sin2^H5 cos2^H, and sin2^H are their values at the resonance position, A{tres) and i?(ti.es)- 
In this case we can replace ^ -|- -B by {5m'^/2Eji) cos 2^ wherever it occurs. Eqs. (20a & 
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b) with = can then be written as: 

tan 2^ 

sin2^H = , (29a) 
V(l - ^/^r)2 + tan^ 2^ 

cos 2^H = ^-E/E^ (296) 

V(l - £;/Er)2 + tan^ 29 

. tan 26* , 

sin2^H = , (29c) 
V(l + ^/£^r)2 + tan2 29 

cos 2^H = ^ ^^^^ . (29d) 

V(l + £^/£^r)2 + tan2 2^ 

(2.) We employ these approximations for the full effective mixing angle to obtain 
estimates for B in Eq. (22a). 

(3.) So far we have not specified the resonance position. We now use A and the 
estimate of B from Step (2) to estimate the resonance position through {5m'^/2E^) cos 29 = 
A + B. Note that A and B are position dependent. 

(4.) With the resonance position from Step (3) we use Eq. (22b) to estimate B^t- 

(5.) With this estimate for B^j. we now can re-estimate the full effective mixing angle 
using 

. {6my2E,)sm29 + B^r , 

sm2^H = — , (30a) 

y/[{5m^/2E,,) - (5mV2£;R)]2cos2 29 + [(Sm^ /2E,,) sin29 + B^r]^ 

y/[{5m^/2E^) - {5m^/2E^)]'^cos^ 29 + [iSm^/2E^)sm29 + B^r]^ 

. {Smy2E,)sm29-Ber . 

sm2^H = — , , (30c) 

^/[{5m^ /2E^) + {Sm^/2En)]^ cos^ 29 + [(Sm^/2E^) sin 2^ - B^r]^ 

y/[{5m^ /2E^) + (5m^/2En)]^ cos^ 29 + [{5w? /2E^) sin 29 - B^r? 
(6.) We iterate by returning to Step (2) and re-evaluating B. 

This procedure must be continued until S, B^t-, ^'h and the resonance position (tres) 
converge. Because of the dependence of B and B^t on the flavor evolution histories of 
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all neutrinos in the ensemble, convergence of this procedure is, in general, problematic. 
However, if neutrino flavor evolution is adiabatic then the complication of prior histories 
is eliminated and the above procedure converges rapidly for the conditions which obtain 
in the region above the neutrino sphere in Type II supernovae. For nonadiabatic neutrino 
flavor evolution the above procedure, though more laborious, still gives good estimates of 
the effects of the neutrino background. We shall begin by discussing the case of adiabatic 
neutrino flavor evolution. 

Ilia.) Adiabatic Neutrino Flavor Evolution 

Consider the flavor evolution of antineutrinos. It is generally true everywhere above 
the neutrino sphere that the contributions of the electrons and neutrinos satisfy A + B > 0. 
This is true because ng is everywhere greater than the net neutrino number densities for 
any neutrino flavor [5,6]. For an antineutrino emitted from the neutrino sphere in the jpg) 
flavor eigenstate, it is evident that |aie(t)P ~ 1 and |aiT-(t)p ft! for all t. The effective 
mass-squared difference for two antineutrino mass eigenstates always increases with density 
and there is no mass level crossing. The adiabatic approximation for the evolution of the 
antineutrino mass eigenstates is always good. 

The situation is more complicated for neutrinos. However, the approximation of 
adiabatic evolution of the neutrino mass eigenstates is a particularly simple case to treat in 
the supernova. A neutrino created in a flavor eigenstate Iva) at the neutrino sphere is very 
nearly in a mass eigenstate because of the large electron number density there. Subsequent 
adiabatic evolution then implies that, for example, |aie(t)P = and |aiT-(t)p = 1 for all 
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t (likewise, |a2e(^)P = 1 and |a2r(i)P = for all t). In this case the expressions for the 
neutrino background contributions, Eqs. (22a & b), become 

B . /(I - oos.,)[<«2.„ - *0 + C0S25„«)(*, - (31a) 

« /(I - cos«,)[si„2«„W(*.. - + .in2«-„(0(*, - </n,J]. (316) 

The evaluation of Eqs. (31a & b) for particular neutrino distribution functions is 
straightforward so long as the adiabatic approximation obtains. To begin with, consider 
the computation of B from Eq. (31a) in the limit where Ser = 0. The result so obtained 
will be valid if we can later show that l-Berl <S {5m? /2E^) sin 2^. 

With the approximation that Ser is small the integrals over the neutrino distribution 
functions (in^^, rfn^^, duy^, and duy^ can be separated into an angular part and an energy 
part. This is due to the fact that when B^.^ is small and essentially become functions 
of energy alone. For a radially propagating neutrino, the angular part of the integral in 
Eq. (31a) then becomes. 



1 /"^° 1 
(1 - cos^q)^ = ~ "^^'^^ "'""^^^ " 4^^ ~ Vl-(i?./r)2]2. 



(32a) 



In this equation r is the radius of the point at which we evaluate B and is the polar angle 
of the limb of the neutrino sphere as seen from this point. Frequently we are interested in 
regions sufficiently distant from the neutrino sphere that we can take r ^ R^. In this limit, 
the radial neutrino path to the point at radius r is a good representation of all neutrino 
paths to that point, and we can approximate 

.rffiq 1 Rf 



f rffiq 1 Rt 

/ (1 - COSb'a)— fti 777 

J ^ ^' Att 16 H 
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It is obvious from this expression that the diagonal contribution of the neutrino background 
is sensitive to position. 

The integration of the remaining energy dependent terms in Eq. (31a) is simple if 
we employ the approximate energy spectra in Eq. (12c). The energy part of Eq. (31a) is 
then, 

J cos 2^h/.„ {E,^ )dE,^ ^ {6, E^/T,^ ) , (33a) 
J cos2enUAEi>JdE^^ « F^{e,En/T^J, (336) 

where we define the neutrino spectral integrals as, 



(34a) 



F^{e,xn) = -i- r ^ + ^/^^ ^! ^3-. (34^) 

' i^2(0)7o V(l + a^/^R)' + tan^ 29 exp(a;) + 1 ^ ^ 

Clearly, for tan 2^ <^ 1, Fi^{6,x^) 1. Here E'r is the energy corresponding to a neutrino 
at resonance at radius r. 

With these definitions, and for small Sgr, we can reduce Eq. (31a) for B to. 



fi _ . /i _ /^2l2 
B ^ -V2G^^- ^ ^ J [<F.(^, E^/T^J - nlF,{0, E^/T,^) 

+ nlF,{e, En/T,^) - n%F,{e, Ej,/T,^)], (35) 



where , , , and are the appropriate neutrino or antineutrino number densities 
at the neutrino sphere as in Eq. (12b). This zero-order expression for B is to be used in 
Step (2) in the iterative procedure outlined above. To proceed further requires that we 
estimate B^r- 
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The angular integration for Eq. (31b) is the same as for Eq. (31a). In performing 
the angular integration in Eq. (31b) we will again assume that is small. The energy 
dependent integrals in Eq. (31b) can be written as: 



(36a) 



j sm2euU^{E^JdE^^ ^ Gi,{e,E^/T^J. (366) 
In like manner to Eqs. (34a & b) we define, 

G.(^,-a)-^y^ /n-./.,)^ + tan-2. exp(.) + l -^-' ^''^^ 



^/{T^^^^J^^^yTt^^^^ exp(a;) + 1 

G (9 X ) = r tan2g x^ 

'^^^ ^2(0)70 V(l + x/xn)'' + tan2 20 exp(a;) + 1 ^' 



(376) 



y(rq:^7^^j2qpt8m22^ exp(a;) + 1 

where the notation is as in Eqs. (34a & b). 

Finally, we can utilize Eqs. (36a-37b) to give an approximate expression for B^t, 

+ nlG^{9,En/T^J-nlG,{9,En/T^^)]. (38) 

The notation in this equation is the same as in Eq. (35). The approximations for B and 
in Eqs. (35) and (38), respectively, are valid when, |Ser|/((^?77,^/2£^R) sin2^ <^ 1. 
Note that the integrand in the expression for Gy{9,x-R) in Eq. (37a) contains a fac- 
tor, sin^H ~ tan 2^/^(1 - x/xr)^ + tan^ 29 = tan 2^/^/(1 - E^, / E^)^ + tan^ 29, which is 
sharply peaked at E^, = for small vacuum mixing angles. In Fig. 5 we plot sin2^H as 
a function of E^/E^ for three values of the vacuum mixing angle. The dotted line in this 
figure corresponds to tan 2^ = 10~^. The dashed line corresponds to tan 2^ = 10~^, while 
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the solid line is for tan 2^ = 0.1. Since the factor sin 2^h appears in the integration over the 
neutrino energy spectrum we can see easily that the smaller the vacuum mixing angle, the 
smaller will be the fraction of the total number density of neutrinos which contribute to 
Ser- The physical interpretation of this is clear: the neutrinos which make the largest con- 
tribution to the off-diagonal neutrino background terms are those which have the largest 
full effective mixing angles at the position under consideration at radius r. These are the 
neutrinos which have energies close to E^. 

With the iterative procedure outlined above we can estimate B and for adiabatic 
neutrino flavor evolution in both the shock re- heating epoch and the hot bubble/r-process 
nucleosynthesis epoch. Refs. [5] and [6] give detailed expositions of the expected neutrino 
emission parameters for these epochs. Typical neutrino luminosities for the shock re- 
heating epoch at tpB ~ 0.15 s (see the discussion in Ref. [5]) are L^^ ^ Lj^^ f» -^^^^c/^) ~ 
Lii^(i^) ~ 5 X 10^^ erg s~-^. The neutrino sphere radius at this epoch is R^^ ^ 50 km, 
while the average neutrino energies are (E,^^) fa 9 MeV, {Ej^J 12 MeV, and (E,^^^^^) 
{Ey^^^^) pa 20 MeV. By contrast, in the later hot bubble/r-process nucleosynthesis epoch 
(tpB fa 5 s) the neutrino liminosities are Lj^^ ^ Ljj^ L^^^^^ Lj^^^^^ lO^-*^ erg s~^, while 
the neutrino sphere is at radius R,^ 10 km. The average neutrino energies for this epoch 
are {E^^) ^ 11 MeV, {E^^J ^ 16 MeV, and {E^^^^^) ^ {E^^^^^) ^ 25 MeV. 

As Ref. [5] shows, for a substantial enhancement in shock re-heating Uj. (or Z/-^) 
neutrinos with energies E,^ Ri 35 MeV must be efficiently transformed into Ve neutrinos in 
the region behind the stalled shock. Ref. [6] shows that neutrinos with 

26 



MeV are the most important in determining the electron fraction, Y^,, in the hot bubble/r- 
process nucleosynthesis epoch. 

For shock re-heating enhancement we must have adiabatic transformation of neutrinos 
with energies E^, Ri 35 MeV. In the hot bubble/r-process nucleosynthesis epoch adiabatic 
transformation is not necessary to drive the material too proton rich for r-process nucle- 
osynthesis to occur (Yg > 0.5). In fact Ref. [6] shows that i^t{h) ^ flavor conversion 
efficiencies as small as ~ 30% for neutrinos with energies E^, 25 MeV will suffice to drive 
Ye > 0.5. Nevertheless, for large enough vacuum mixing angles, adiabatic transformation 
of neutrinos with Ei^ 25 MeV will occur in some regions of the {Sm'^, sin^ 29) plot (Fig. 
2 in Ref. [6]). 

Consider adiabatic neutrino flavor conversion specifically for Ei, = 35 MeV in the 
shock re- heating epoch and Ei^ = 25 MeV in the hot bubble/r-process nucleosynthesis 
epoch. For comparison, we first present values of 5m^ and sin^ 29 which give an adiabaticity 
parameter 7 = 3 for the bare electron number density distributions relevant for these 
epochs. In Figs. 6 and 7, the solid contour lines for 7 = 3 correspond to these values of 
Sm^ and sin^ 29 for the representative conditions in the shock re- heating and hot bubble/r- 
process nucleosynthesis epochs, respectively. 

It should be noted that the adiabatic approximation will be valid over the whole 
range of neutrino energies implicit in the neutrino distribution functions entering into the 
expressions for B and Ser- Neutrinos with energies Ei, < E^_ will propagate through 
resonances prior to reaching the resonant position for the specific example neutrino energy 
under discussion (either £^r, = 35 MeV or E^ = 25 MeV). It is a general feature of the 
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density scale height above the neutrino sphere that the neutrinos with energies < 
will experience adiabatic flavor evolution through their resonances as long as neutrinos 
with energy E'r go through the resonance adiabatically [5, 6]. Background neutrinos with 
energies E^ > £'r will not have gone through resonances and therefore evolve adiabatically 
prior to arriving at the resonance position for a neutrino with energy E^^. We conclude 
that the expressions for B and Ber in Eqs. (31a & b) are appropriate for the example 
under consideration. 

Using the iterative procedure outlined above we can calculate the true adiabatic pa- 
rameter, 7(tres): including the neutrino background contributions. We show the new con- 
tour lines for 7 = 3 as dotted lines in Figs. 6 and 7 for the respective epochs. We can 
easily see that the neutrino background has a completely negligible effect on adiabaticity 
at resonance along the solid 7 = 3 contour line in Fig. 6. The new contour line for 7 = 3 
in Fig. 6 is indistinguishable from the one calculated for the bare electron number density. 
The new contour line for 7 = 3 in Fig. 7 moves a little bit to the right of the solid line, 
but the neutrino background effects are also evidently small. 

Any neutrino mixing parameters Sm? and sin^ 29 which are to the right of the 7 = 3 
contour lines in Figs. 6 and 7 correspond to larger values of 7 for the specific example 
neutrino energies under discussion. For a given Sm'^ the ratio \Ber\/ (Sm'^ /2En) sin 2^ 
will decrease as sin^ 29 and, hence, 7 increases. The off-diagonal neutrino background 
contribution will have a negligible effect on neutrino flavor conversion everywhere to the 
right of the contour lines in Figs. 6 and 7. Likewise, B is roughly constant for a given Sm^ 
as sin^ 29 and, hence, 7 is increased. The diagonal contribution of the neutrino background 
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produces a negligible alteration in the computed flavor conversion efficiencies everywhere 
to the right of the contour lines in Figs. 6 and 7. 

We have also examined adiabatic neutrino flavor conversion in supernovae for a range 
of neutrino energies. We can conclude that the neutrino background, specifically B and 
Ser, will not result in any modification of the results of Refs. [5] and [6] whenever adiabatic 
neutrino flavor evolution is at issue. 

Illb.) Nonadiabatic Neutrino Flavor Evolution 

The effects of the neutrino background on nonadiabatic neutrino flavor evolution in 
the region above the neutrino sphere are potentially more signiflcant than are the neutrino 
background effects on adiabatic neutrino flavor evolution. In general, the evaluation of 
B and Bgj. from Eqs. (22a & b) is considerably more complicated when neutrino flavor 
evolution is nonadiabatic than it is when the adiabatic limit for neutrino flavor evolution 
obtains. 

A neutrino of energy En, nonadiabatically going through a resonance at a point above 
the neutrino sphere, experiences a neutrino background effect which depends on the prior 
histories of all the neutrinos in the ensemble which are passing through the resonance 
region. In this case, we cannot argue that background neutrinos with E^, < E^ go through 
resonances adiabatically. The flavor evolution for background neutrinos with E,^ > E^i 
can still be considered adiabatic for the purposes of calculating B and Ser, since these 
neutrinos will not yet have gone through resonances when they are in the resonance region 
for energy Er. 
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In Fig. 8 we graphically illustrate the difficulties inherent in computing B and Bg^- 
from Eqs. (22a & b) for nonadiabatic neutrino flavor evolution. In this figure we show 
the radial path of a neutrino with energy E'r. The resonance position for this neutrino 
is the point labeled RES(£'r). The path for a neutrino of energy Eb representative of 
the neutrino background at the point RES(£^r) is labeled by Eb- If Eb < -Er, then the 
neutrino on the path labeled by Eb presumably propagated through a resonance of its own 
prior to reaching position RES(£?r.). The resonace position for the background neutrino is 
labeled RES(£^b)- Whether or not this background neutrino experiences flavor conversion 
at RES(£'b) depends, in turn, on the flavor evolution histories of the background neutrinos 
which pass through this point. The paths for some of these "secondary" background 
neutrinos are shown in Fig. 8. 

As we can see from Fig. 8, an exact calculation of the neutrino background contribu- 
tions requires us to simultaneously follow the flavor evolution histories of neutrinos with 
different energies on all possible neutrino paths above the neutrino sphere. This could be 
done in a Monte Carlo calculation. However, there is a simpler alternative if we make 
note of the following two facts. First, we are most interested in regions which are far away 
from the neutrino sphere. The region for r-process nucleosynthesis in the hot bubble is 
located at radii r > 4:Rj^. So the polar angles for neutrino paths to a point in this region 
lie in a narrow range around = 0. In addition, at a point close to the neutrino sphere 
where the polar angles for the relevant neutrino paths can be significantly different from 
zero, the electron number density is so high that neutrino background effects can be safely 
ignored. Therefore, we can make an approximation and take the flavor evolution history of 
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a radially propagating neutrino (^q = 0) as representative of the flavor evolution histories 
of all neutrinos with the same energy. 

The flavor evolution history of radially propagating neutrinos for a given set of Sm^ 



{!'.) We numerically represent the neutrino energy spectrum with a grid of energy 
bins. These energy bins cover a neutrino energy range of 1-100 MeV. Typically our nu- 
merical calculations employ ~ 200 energy bins. Since neutrinos with lower energies go 
through resonances first, we start the calculations at the lower end of the energy grid. 

(2'.) For the particular grid point (neutrino energy bin) at neutrino energy Ei,, we 
use the iterative procedure outlined at the beginning of this section to locate the resonance 
position, tj-ea{E^), for this particular neutrino energy Ei^. As a byproduct of this iterative 
procedure, we will obtain the corresponding neutrino background contributions B and Bgj. 
at this position trea{Ei^). The evaluation of B and Ber in this case is quite similar to that 
for the case of adiabatic neutrino fiavor evolution, except that here we must use Eqs. (22a 
& b) together with the flavor evolution histories of neutrinos with energies lower than E,^. 

(3'.) Using the resonance position, tres{Eu), and the corresponding neutrino back- 
ground contributions B and Bgr from Step (2'), we can evaluate the Landau-Zener proba- 
bility Plz(-Ei/) (Eq. [28]) for a neutrino with energy E,^ to jump from one mass eigenstate 
to the other in the course of transversing the resonance region. 

(4'.) The flavor evolution history of i^g neutrinos with energy Ejy is then approximated 

as 



and sin 



^ 29 can then be calculated with the following procedure: 





otherwise. 



(39) 
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Likewise, the flavor evolution history of neutrinos with energy E^, is approximated as 



In the above two equations, the evolutionary parameter t increases away from the neutrino 
sphere. These approximations for the neutrino flavor evolution history, together with Eqs. 
(22a & b), are then used in the iterative procedure in Step (2') to locate the resonance 
position and calculate the corresponding neutrino background contributions for neutrinos 
with energies higher than E^^. 

At the end of the above procedure, we will have obtained the approximate flavor 
evolution histories for all the neutrino energies on the energy grid. This information then 
can be used to calculate the electron fraction Yg in the r-process nucleosynthesis region as 
described in Ref. [6]. We present the new Ye = 0.5 line, including the neutrino background 
effects, as a dotted contour line on the {Sw? , sin^ 29) plot in Fig. 9. The original = 0.5 
line in Fig. 2 of Ref. [6] is shown as the solid contour line in Fig. 9. To the right of the 
Yg = 0.5 line, the material will be driven too proton rich for r-process nucleosynthesis to 
occur in the hot bubble. 

By examining the two contour lines in Fig. 9, we can draw two conclusions. First, with 
a proper treatment of the neutrino background effects, we see that r-process nucleosynthe- 
sis in the hot bubble remains a sensitive probe of the flavor-mixing properties of neutrinos 
with cosmologically signiflcant masses. In fact, inclusion of the neutrino background con- 
tributions results in a small modiflcation of the original Yg = 0.5 line for dm? = 4 eV^ to 
dm^ — 10"^ eV^. Furthermore, after we take into account the neutrino background con- 
tributions, it is evident that the range of neutrino vacuum mass-squared difference bvrP' 




(40) 
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probed by r-process nucleosynthesis is extended down to Sm? < 2 eV^. The reason for this 
extension can be found in the nonhnear nature of neutrino flavor transformation in the 
presence of a neutrino background. 

Close to the neutrino sphere where little neutrino flavor transformation has occurred, 
the number density of Ve neutrinos is larger than that of Vr neutrinos. This is because 
the luminosities for i/^ and i/r are approximately the same, but the average i/r neutrino 
energy is much higher (cf. Eq. [12b]). However, with neutrino flavor transformation, more 
Ve neutrinos are transformed into Vr neutrinos than Vr neutrinos are transformed into 
neutrinos. This is because there are more low energy neutrinos and only low energy 
neutrinos are very efficiently transformed for the parameters along the dotted contour line 
in Fig. 9. Because of the nonlinear evolution of the neutrino background, the diagonal 
contribution B evolves from a positive value for positions close to the neutrino sphere to 
a negative value for positions far away from the neutrino sphere. Neutrinos with Srn^ < 2 
eV^ and energies over a broad range will tend to have resonances far enough out that the 
diagonal contributions will satisfy B < 0. For a given Sm^ and a given energy Ei^, the 
resonance position will lie closer to the neutrino sphere for the case B < than it would 
for the case where no neutrino background is present (cf. Eq. [26]). 

As Ref. [6] discusses. Ye and, hence, r-process nucleosynthesis are sensitive to neutrino 
flavor conversion only when resonances occur inside the weak freeze-out radius. The weak 
freeze-out radius is the radius beyond which typical i^e and Pg capture rates are small 
compared to the material expansion rate. When B <0, the resonances for given 5m? are 
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drawn in toward the neutrino sphere. Hence, we find that the Yg = 0.5 fine drops to lower 
values of Sm^ in the presence of a neutrino background. 



IV. Conclusions 

We have calculated neutrino flavor transformation in the region above the neutrino 
sphere in Type II supernovae including all contributions from the neutrino background. In 
particular, we have examined the neutrino background effects on both cases of adiabatic 
and nonadiabatic neutrino flavor evolution. In the case of adiabatic neutrino flavor evo- 
lution, which is most relevant for supernova shock re-heating, we find that the neutrino 
background has a completely negligible effect on the range of vacuum mass-squared dif- 
ference, 5m^, and vacuum mixing angle, 9, or equivalently sin^ 2^, required for enhanced 
shock heating. In the case of nonadiabatic neutrino fiavor evolution relevant for r-process 
nucleosynthesis in the hot bubble, we find that r-process nucleosynthesis from neutrino- 
heated supernova ejecta remains a sensitive probe of the mixing between a light and a 
i^t{h) with a cosmologically significant mass. The modification of the (5m^, sin^ 26*) pa- 
rameter region probed by r-process nucleosynthesis due to the neutrino background effects 
is generally small. The nonlinear nature of neutrino fiavor transformation in the presence 
of a neutrino background actually extends the sensitivity of r-process nucleosynthesis to 
smaller values of 5m^. 

In general, we find that a proper account of neutrino background effects leads to 
no modification in the overall qualitative conclusions of Refs. [5] and [6]. At the early 
epochs of the post-core-bounce supernova environment (tpB < 1 s), we find that the 
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characteristically large electron number densities and large density scale heights determine 
the phenomenon of neutrino flavor transformation. Even at the later epochs associated 
with r-process nucleosynthesis, the effects of the neutrino background on neutrino flavor 
evolution are small. 
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Figure Captions 



Fig. 1 A generic Feynman graph for v^-e scattering. 

Fig. 2 The geometrical arrangement of a neutrino sphere with radius Rv^ a point above 
the neutrino sphere at radius r, and various neutrino paths. 

Fig. 3 Generic Feynman graphs for neutrino-neutrino exchange-scattering processes. Fig. 
3a is for Vf.-Vf. scattering and Fig. 3b is for Vt-Vt scattering. 

Fig. 4 Graphic representation for off-diagonal contributions from the neutrino back- 
ground. 

Fig. 5 The zero-order expression for sin2^H as a function of Ej^/Eyi for three different 
vacuum mixing angles. The dotted line corresponds to tan 2^ = 10""^. The dashed line 
corresponds to tan 2^ = 10~^, while the solid line is for tan 2^ = 0.1. 
Fig. 6 Contour lines for 7 = 3 on the (5m^, sin^ 29) plot for the shock re-heating epoch. 
The solid contour line is calculated for the bare electron number density. The dotted 
line, which cannot be distinguished from the solid line in this case, is calculated with the 
neutrino background contributions. 

Fig. 7 As in Fig. 6, but for the hot bubble/r-process nucleosynthesis epoch. 
Fig. 8 Illustration of the difficulties inherent in computing the neutrino background 
contributions B and -Ber for the case of nonadiabatic neutrino flavor evolution. The radial 
path of a neutrino with energy and resonance position RES(£^r,) is shown. The path for 
a neutrino of energy E^ representative of the neutrino background at position RES(£^r) 
is shown together with its resonance position RES(i?B)- Paths for background neutrinos 
at position RES(£?b) are also shown. 
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Fig. 9 Contour lines for = 0.5 are shown on the {5m?, sin^ 20) plot. The solid line is 
the same as the Yg — 0.5 line in Fig. 2 of Ref. [6], whereas the dotted line is calculated 
with the full neutrino background contributions. 
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